library(Seurat)
library(dplyr)
library(tibble)
library(cowplot)
Data loading and processing

General exploratory analysis


U/L branch marker genes
Plot the U and L branch marker genes.
#DE genes between T1T2T3 and T4T5 in the 10x-180831 data.
markers_u_l <- read.table('output/markergenes/180831/markers_10x-180831_upperbranch_lowerbranch_negbinom', sep='\t', header=T)
cannot open file 'output/markergenes/180831/markers_10x-180831_upperbranch_lowerbranch_negbinom': No such file or directoryError in file(file, "rt") : cannot open the connection
How do these genes look in the 180831 data?

And how are they expressed in the Vijay data?

plots <- FeaturePlot(vijay, features=c(as.vector(markers_u$gene)[1:10], as.vector(markers_l$gene)[1:10]), pt.size=1, combine=F)
The following requested variables were not found: RP11-572C15.6
plot_grid(plotlist=plots, ncol=2)

Some U branch markers are not expressed at all (SCD, PLIN4, GDP1). The other U branch markers are expressed in P2 (FABP4), in I1 and I3 and E1-E3.
Seurat integration with 10x-180831 data
Comparison with 10x-180831 data
check marker genes from their paper in 10x-180831 data
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3IgbWVzc2FnZT1GfQpsaWJyYXJ5KFNldXJhdCkKbGlicmFyeShkcGx5cikKbGlicmFyeSh0aWJibGUpCmxpYnJhcnkoY293cGxvdCkKYGBgCgpEYXRhIGxvYWRpbmcgYW5kIHByb2Nlc3NpbmcKCmBgYHtyfQpmaWxlLmRhdGEgPC0gIi9uZnNkYXRhL3Byb2plY3RzL3RpbXNoZWwvc2Mtc2NoZWVsZV9sYWJfYWRpcG9zZV9mbHVpZGlnbV9jMS9kYXRhLXZpamF5L0dTRTEyOTM2M19TVkZfTm9ybWFsaXNlZF9EYXRhLnR4dC5neiIKZGYuZGF0YSA8LSByZWFkX3RzdihmaWxlLmRhdGEpICU+JSByZW5hbWUoZ2VuZT1YMSkjIFdPUktTIGFuZCBwYXJzZWQgY29ycmVjdGx5LiBXYXJuaW5nIG1lc3NhZ2U6IE1pc3NpbmcgY29sdW1uIG5hbWVzIGZpbGxlZCBpbjogJ1gxJyBbMV0gCmRmLmRhdGEgPC0gZGF0YS50YWJsZTo6ZnJlYWQoZmlsZS5kYXRhLCBuVGhyZWFkPTI0LCBzaG93UHJvZ3Jlc3M9VCkKZGYuZGF0YSA8LSBkZi5kYXRhICU+JSBjb2x1bW5fdG9fcm93bmFtZXModmFyPSJWMSIpCgojIyMgTWV0YS1kYXRhCmZpbGUubWV0YWRhdGEgPC0gIi9uZnNkYXRhL3Byb2plY3RzL3RpbXNoZWwvc2Mtc2NoZWVsZV9sYWJfYWRpcG9zZV9mbHVpZGlnbV9jMS9kYXRhLXZpamF5L0dTRTEyOTM2M19jZWxsX21ldGFkYXRhLnR4dC5neiIKZGYubWV0YWRhdGEgPC0gcmVhZF90c3YoZmlsZS5tZXRhZGF0YSkKZGYubWV0YWRhdGEgPC0gZGYubWV0YWRhdGEgJT4lIHJlbmFtZShjZWxsX2JhcmNvZGU9WDEsIGRpYWJldGVzPUNvbmRpdGlvbiwgZGVwb3Q9VGlzc3VlLCBjZWxsdHlwZT1DbHNJRCkKIyAqVE9ETyo6IHJlbmFtZSBJQS0tPnZpc2M7IFNDLS0+c3VicQpkZi5tZXRhZGF0YSA8LSBkZi5tZXRhZGF0YSAlPiUgY29sdW1uX3RvX3Jvd25hbWVzKHZhcj0iY2VsbF9iYXJjb2RlIikKCiMjIyBDcmVhdGUgb2JqZWN0CnZpamF5IDwtIENyZWF0ZVNldXJhdE9iamVjdChjb3VudHM9ZGYuZGF0YSwgbWV0YS5kYXRhPWRmLm1ldGFkYXRhLCBwcm9qZWN0PSJ2aWpheSIpICMgV2FybmluZzogRmVhdHVyZSBuYW1lcyBjYW5ub3QgaGF2ZSB1bmRlcnNjb3JlcyAoJ18nKSwgcmVwbGFjaW5nIHdpdGggZGFzaGVzICgnLScpCnZpamF5W1siUk5BIl1dQGRhdGEgPC0gdmlqYXlbWyJSTkEiXV1AY291bnRzICMgY29weSBzbG90IGJlY2F1c2UgZGF0YSBpcyBhbHJlYWR5IG5vcm1hbGl6ZWQKIyBOQjogZG8gbm90IHBlcmZvcm0gZnVydGhlciAobG9nKW5vcm1hbGl6YXRpb24gYXMgdGhpcyBoYXMgYWxyZWFkeSBiZWVuIGRvbmUKCnZpamF5IDwtIFNjYWxlRGF0YSh2aWpheSkKdmlqYXkgPC0gRmluZFZhcmlhYmxlRmVhdHVyZXModmlqYXkpCnZpamF5IDwtIFJ1blBDQSh2aWpheSkKRWxib3dQbG90KHZpamF5LCByZWR1Y3Rpb249J3BjYScpCmBgYAoKYGBge3J9CnZpamF5IDwtIFJ1blRTTkUodmlqYXksIGRpbXM9MToxMSkKc2F2ZVJEUyh2aWpheSwgJ291dHB1dC9zZXVyYXRfb2JqZWN0cy92aWpheS92aWpheS5yZHMnKQpgYGAKCiNHZW5lcmFsIGV4cGxvcmF0b3J5IGFuYWx5c2lzCgpgYGB7ciBmaWcuaGVpZ2h0ID0gMTAsIGZpZy53aWR0aCA9IDEyLCBmaWcuYWxpZ24gPSAiY2VudGVyIn0KI3Bsb3QgY2VsbCB0eXBlcyArIG92ZXJsYXllZCB3aXRoIG1ldGEgZGF0YQpwbG90X2dyaWQoCiBUU05FUGxvdCh2aWpheSwgZ3JvdXAuYnk9J2NlbGx0eXBlJywgbGFiZWw9VCksIAogVFNORVBsb3QodmlqYXksIGdyb3VwLmJ5PSdkZXBvdCcpLAogVFNORVBsb3QodmlqYXksIGdyb3VwLmJ5PSdkaWFiZXRlcycpCikKYGBgCgpgYGB7ciBmaWcuaGVpZ2h0ID0gNSwgZmlnLndpZHRoID0gMTIsIGZpZy5hbGlnbiA9ICJjZW50ZXIifQojcGxvdCBjZWxsIHR5cGVzICsgb3ZlcmxheWVkIHdpdGggbWV0YSBkYXRhCnBsb3RfZ3JpZCgKIEZlYXR1cmVQbG90KHZpamF5LCBmZWF0dXJlcz0nbkZlYXR1cmVfUk5BJyksIAogRmVhdHVyZVBsb3QodmlqYXksIGZlYXR1cmVzPSduQ291bnRfUk5BJykKKQpgYGAKCiNVL0wgYnJhbmNoIG1hcmtlciBnZW5lcwoKUGxvdCB0aGUgVSBhbmQgTCBicmFuY2ggbWFya2VyIGdlbmVzLiAKCmBgYHtyfQpkYXRhXzE4MDgzMSA8LSByZWFkUkRTKCdvdXRwdXQvc2V1cmF0X29iamVjdHMvMTgwODMxLzEweC0xODA4MzEtUzMnKQojREUgZ2VuZXMgYmV0d2VlbiBUMVQyVDMgYW5kIFQ0VDUgaW4gdGhlIDEweC0xODA4MzEgZGF0YS4KbWFya2Vyc191X2wgPC0gcmVhZC50YWJsZSgnb3V0cHV0L21hcmtlcmdlbmVzLzE4MDgzMS9tYXJrZXJzXzEweC0xODA4MzFfdXBwZXJicmFuY2hfbG93ZXJicmFuY2hfbmVnYmlub20nLCBzZXA9J1x0JywgaGVhZGVyPVQpCm1hcmtlcnNfdSA8LSBtYXJrZXJzX3VfbFtvcmRlcigtbWFya2Vyc191X2wkYXZnX2xvZ0ZDKSxdCm1hcmtlcnNfbCA8LSBtYXJrZXJzX3VfbFtvcmRlcihtYXJrZXJzX3VfbCRhdmdfbG9nRkMpLF0KYGBgCgpIb3cgZG8gdGhlc2UgZ2VuZXMgbG9vayBpbiB0aGUgMTgwODMxIGRhdGE/CgpgYGB7ciwgZmlnLmhlaWdodCA9IDEyLCBmaWcud2lkdGggPSAxMiwgZmlnLmFsaWduID0gImNlbnRlciJ9CnBsb3RzIDwtIEZlYXR1cmVQbG90KGRhdGFfMTgwODMxLCBmZWF0dXJlcz1jKGFzLnZlY3RvcihtYXJrZXJzX3UkZ2VuZSlbMToxMF0sIGFzLnZlY3RvcihtYXJrZXJzX2wkZ2VuZSlbMToxMF0pLCBwdC5zaXplPTEsIGNvbWJpbmU9RikKcGxvdF9ncmlkKHBsb3RsaXN0PXBsb3RzLCBuY29sPTQpCmBgYAoKQW5kIGhvdyBhcmUgdGhleSBleHByZXNzZWQgaW4gdGhlIFZpamF5IGRhdGE/CgpgYGB7cn0KVFNORVBsb3QodmlqYXksIGdyb3VwLmJ5PSdjZWxsdHlwZScsIGxhYmVsPVQpCmBgYAoKCmBgYHtyLCBmaWcuaGVpZ2h0ID0gNTAsIGZpZy53aWR0aCA9IDEyLCBmaWcuYWxpZ24gPSAiY2VudGVyIn0KcGxvdHMgPC0gRmVhdHVyZVBsb3QodmlqYXksIGZlYXR1cmVzPWMoYXMudmVjdG9yKG1hcmtlcnNfdSRnZW5lKVsxOjEwXSwgYXMudmVjdG9yKG1hcmtlcnNfbCRnZW5lKVsxOjEwXSksIHB0LnNpemU9MSwgY29tYmluZT1GKQpwbG90X2dyaWQocGxvdGxpc3Q9cGxvdHMsIG5jb2w9MikKYGBgCgpTb21lIFUgYnJhbmNoIG1hcmtlcnMgYXJlIG5vdCBleHByZXNzZWQgYXQgYWxsIChTQ0QsIFBMSU40LCBHRFAxKS4gVGhlIG90aGVyIFUgYnJhbmNoIG1hcmtlcnMgYXJlIGV4cHJlc3NlZCBpbiBQMiAoRkFCUDQpLCBpbiBJMSBhbmQgSTMgYW5kIEUxLUUzLiAKCiNTZXVyYXQgaW50ZWdyYXRpb24gd2l0aCAxMHgtMTgwODMxIGRhdGEKCmBgYHtyfQoKYGBgCgoKI0NvbXBhcmlzb24gd2l0aCAxMHgtMTgwODMxIGRhdGEKCmBgYHtyfQpjaGVjayBtYXJrZXIgZ2VuZXMgZnJvbSB0aGVpciBwYXBlciBpbiAxMHgtMTgwODMxIGRhdGEKYGBgCgo=